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Abstract 

We study the interplay between nematic order and superconductivity, motivated by a recent 
experiment on FeSe observing strongly distorted vortex shapes (Song et al, Science 332, 1410 
(2011)). We show that the nematic order strongly enhances the anisotropy in the superconducting 
coherence length, beyond that expected from considerations of the Ginzburg-Landau theory. We 
obtain universal functions describing the coupling between the nematic order and superconductiv- 
ity, and discuss connections of our results to the experiments. 



I. INTRODUCTION 



The presence of multiple order parameters is one of the most interesting features of 
strongly correlated systems. For example, a number of order parameters appear in both 
the copper and iron based high temperature superconductors, including superconductivity 
and spin density wave order (SDW). Among them, the electronic nematic order parameter, 
measuring spontaneous symmetry breaking of lattice rotation symmetry via electron correla- 
tions, has been of particular interest in a series of remarkable experiments)^ In the cuprates, 
nematic order is signaled by the anisotropy of resistivity and bond ordered density^'. In 
the pnictides, orthorhombic structure distortion is ubiquitous, and is confirmed by various 
tools such as neutron scattering^, STJVp^, transportP^Ei', and ARPES 14 . A natural question 
associated with nematic ordering is its interaction with other order parameters. In super- 
conducting materials, we are interested in whether the SC and nematic orders compete or 
attract. In the pnictides, and especially in Bai[Fe2- x Co x ]As2 materials, it was shown that 
the nematic order parameter and SC compete with each other, by observing the suppression 
of structural anisotropy along the onset of SC, the so-called "back-bending"^^ of the phase 
boundary. 

Recently, Song et a/P reported another interesting experiment in the FeSe pnictides ma- 
terial showing striking interplay between nematic order and SC. By observing scanning 
tunneling microscopy (STM) data of the zero bias resonance peak of the electronic density 
of states, Song et al. found that shape of vortex cores was strongly distorted. Even though 
the material has a orthorhombic structure, its structural anisotropy is much smaller than 
the anisotropy of vortices: the reported structural anisotropy is 0.5 percent of the lattice 
constants, but the vortex shape anisotropy is order unity. Song et al. suggested other candi- 
dates, such as orbital ordering, to explain the large vortex anisotropy. Also, the experiments 
showed that the anisotropy is suppressed in strongly SC regions by observing essentially 
isotropic gap functions far from vortices. Thus, it is clear that nematic order and SC com- 
pete each other. Near the vortex core, SC order parameter is suppressed, so nematic order 
is enhanced distorting vortex shapes significantly. 

This paper will present a general theory of competition between the SC and nematic 
orders, as concomitant instabilities of an underlying Fermi liquid. We will not address mi- 
croscopic questions addressed in other literatures^^!, for example, role of orbital physics 

2 



in nematic ordering 24 , but focus on general results on the interplay between the two orders 
that follow from the symmetries of the order parameter and the presence of a Fermi sur- 
face. We will show that the Fermi surface induces strong competition, which is especially 
important in the vortex core, and which can greatly enhance the anisotropy in the vortex 
shape. Non-analyticities associated with the Fermi surface will be shown to induce a large 
coupling between the spatial gradients of the SC order and the nematic order parameter. 
The influence of such a coupling on the vortex shape was studied recently^ in the context 
of a Ginzburg-Landau theory, and our results here offer a rationale for its enhancement. 

Before introducing electrons and their Fermi surface, let us consider a simpler version of 
competing order parameter theory, so-called Ginzburg-Landau theory. 

+ \\^ ( J ) 2 + g GL ^(-d 2 x + d 2 y )i; 

+ ■■■ (1) 

Usual phenomenological constants for the GL theory (r, u, d, y, A, gci) are introduced.^ The 
first and second lines are for the nematic and superconducting order parameters. The third 
line describes the coupling between the two order parameters. The coupling constant A 
characterizes how the two order parameters interact with each other : positive (negative) 
sign means competition (attraction). One of the easiest ways to see the interplay is to observe 
the shift of the nematic order quantum critical point (QCP) shift due to the presence of 
superconductivity: 

Sc (|A|) = Sc (0)-A|A| 2 . (2) 

Here, it is clear that the sign of A determines properties of interplay physics. Its pos- 
itive (negative) value indicates critical point shift toward(away from) the ordered phase, 
which shows competition(attraction) of the two order parameters] 15 * 26 Let us consider an- 
other limit, where the nematic ordering appears with the SC. The SC coherence length 
becomes anisotropic between the x and y directions 

C = \s\- 1/2 (l - g GL 4>/2) Cy = |*r 1/2 (l + W>/2) (3) 



3 



So the gci term describes the anisotropy^ — £ y )/(£ x + £,y) of the SC coherence lengths 
and it is proportial to magnitude of the nematic order linearly. A small nematic order 
parameter implies small difference of coherence lengths in the GL theory. Thus in this 
theory the structural orthorhombicity and nematic ordering cannot explain the strongly 
distorted vortex shape observed by Song et al. 

We will discuss here the modifications of the GL theory described above due to the 
presence of electrons and a Fermi surface. We will show below that the theory with elec- 
trons determines sign of the phenomenological constant A, which turns out to describe the 
competition between the two order parameters. Also, we show that the coupling gcL cor- 
responding to anisotropic coherence length is significantly enhanced by the Fermi surface. 
Quite generally, such effects can be described by the one-particle irreducible vertex function 
in momentum space between two SC and one nematic order, 

(P\ pl-P 2 y ( ) 

, where p XtV ,p,A are for momentum components of SC order parameter, magnitude of the 
momentum, and SC gap magnitude. In the present terms, the GL theory clearly has 

-pGL 2 F GL 2 2 (r\ 

as follows directly from Eq. ([!]). Our main result of this paper is that the vertex function of 
the theory with electrons has a very different behavior 

_ p 2 — p 2 

F^*^ -> const. , V^*^ ~ x ^ 2 y . (6) 

in the p/|A| — > oo limit, which is relevant the near the vortex core region due to SC 
suppression. By comparing the two T^*^ functions at small momentum limit, it is manifest 
that anisotropy is much stronger once we consider electron degrees of freedom. 

The non-analytic term above looks somewhat exotic, but it can be easily understood by 
extensions of the BCS theory in the presence of multiple order parameters. One way to 
understand this is to consider pairing susceptibility with infinitesimal SC pairing. Then, 
pairing susceptibility shows the well known "BCS" logarithmic behavior. 1 27 1 28 1 

Xpairing = (^(p)^(p)) ~ Mp 2 ), (7) 



i]}\x) = cl(x)cUx) is the pairing operator of electrons c a . In the presence of nematic order 
4>, it is clear that there can be a correction to the pairing susceptibility of the form 

Xpainng(0) = (?) V> ip) ) <f> ~ log(p 2 + 4>{vl ~ P?)) (8) 

The vertex function is readily obtained by differentiation with the nematic order parameter, 

l W*4> — 2 — > l 9 J 

which is the result above. Below, we will calculate the full functional dependence of the 
vertex function with two parameters, the momentum and gap magnitudes. 

The main implication of the above considerations is that the linear dependence of the SC 
coherence length on the nematic order obtained in the GL theory is not correct. Even small 
nematic order parameter can show strong enhancement due to its non-analytic influence. 
Thus, we argue that such enhanced competition between SC and nematic order in the 
anisotropic coherence length channel contributes to experimentally observed vortex shape. 

The structure of this paper is as follows. In Sec. [XTJ we set up our model Hamiltonian with 
fermions. We describe how the nematic and the superconducting order parameters couple 
to fermions and introduce our strategy to study the competition effect. For completeness, 
we reproduce the quantum critical point shift of the nematic order parameter with fermions 



under superconductivity in Sec. Ill It is shown that the shift is not significant, in the sense 
that it is of the same order as predicted by the GL theory. Thus, for the QCP shift only, 
the fermions do not play an important role, and the GL theory is enough to describe the 
competition physics. But we show that fermions do play an important role in the anisotropic 



coherence length of the SC order parameter in Sec. |IV| In Sec. [V| we discuss implication of 
our results. 



II. MODEL HAMILTONIAN AND STRATEGY 

In two spatial dimensions, Ising nematic order on the square lattice breaks the four-fold 
rotational symmetry down to two-fold rotation. The Ising nematic order, could be realized 
microcopically by many routes, such as an orbital ordering which may be the case in the 
pnictides. Instead of using orbital ordering however, we focus on a here simpler way to study 
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the nematic order parameter- anisotropic hopping. 

Let us build a model Hamiltonian step by step. A non-interacting Hamiltonian on a 
square lattice is 



where respects the fourfold rotational symmetry. Because a nematic order parameter 
is coupled to an anisotropic hopping, the Hamiltonian with the nematic order parameter 
becomes 



The second term describes an anisotropic hopping, and by condensing the nematic order 
parameter, <p, the four-fold rotational symmetry(or x, y exchange symmetry) is broken down 
to two-fold. In momentum space, it is nothing but order parameter will distort the Fermi 
surface. The third term, H nem , represents the nematic order parameter dynamics. 

If we had used an orbital ordering model, the Fermi surface would contain orbital infor- 
mation, so that each momentum point is tied to the orbital direction. This would modify 
matrix elements on the nematic order parameter in a momentum-dependent way, but not 
crucially modify the Fermi surface physics we are interested in here. 

Next, let us consider the superconducting order parameter. In this paper, we mainly 
focus on the s-wave paired pnictides materials. It is straightforward to generalize our results 
to other cases like ci-wave pairing. After a Fourier transform, the total Hamiltonian is 



where the spin index is suppressed. The H$c is the effective Hamiltonian for the supercon- 
ducting order parameter. Note that, in the continuum limit, the nematic order parameter 
coupling becomes 




(10) 




(11) 




(12) 



4>(cos(k x a) — cos(k y a))c* (k)c(k) — > j(k) <f>c*(k)c(k) 



l(k) = \a\kl-kl). 



(13) 
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The total action is 

Stat = J + A/, + c\(d T + e(k))c k + (j)^{k)c\c k + (V> P 4+ P c -fc + /i.e.) (14) 

The first two terms are for the order parameter dynamics, and one can understand these 
terms as usual (p 4 field theories with one and two components. The remaining terms describe 
the fermion spectrum and its coupling to the two order parameters. Here we absorb coupling 
constants of the Yukawa-type vertices into the order parameters. As expected, the two order 
parameters are coupled to each other by mediating electrons. Below, we study the total 
Lagrangian and how this theory describes the competition between the nematic order and 
the superconductivity. 



To see this, let us introduce some notations. Schematically, the effective action is ex- 
pressed as follows. 



Z = ti(e- H / T ) = [ e-^ff 

J<h,tb 



r e// = J r^Tp*^ + i>^*0 + r H> ^(j>'ip*ijj h — . (15) 

In the previous works^^, the competition between SC and other order parameters was 



described by focusing on terms such as the first two terms in Eqn. 15 There, it was 
concluded that the SC and nematic order do not compete significantly, just as in GL theory, 
compared the competition between SC and spin density wave order (SDW). In this paper, we 
focus on the third term(r^*^) of the above equation, which basically describes anisotropy 
of coherence lengths. 



In this setup, one might argue that one can forget about the fermions in the low energy 
theory, expecting that the SC order parameter induces an energy gap near the Fermi surface. 
Then it seems the Ginzburg-Landau theory of SC and nematic order, Eqn. [TJ would be 
enough to describe the system. However, it turns out that such expectation is too naive 
because we are focusing on physics at non-zero momentum. So we need to keep the fermions 
to see the interplay physics between the order parameters, as shown below. 
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III. NEMATIC QUANTUM CRITICAL POINT SHIFT 



In this section, we consider the nematic QCP shift under the SC. The analysis is almost 
the same as in the previous workp^, and the main difference is that we are consider whole 
Fermi surface, instead of hot-spots. 

Let us consider the quantum critical point as that associated with the onset of nematic 
order, 0, in a metal. This critical point will appear at a value s° of some tuning parameter 
s. We define s so that s < s° is the nematic ordered phase with ^ 0.(see the GL theory- 
above). The value of s° depends on material's microscopic detailes obviously. Then, let 
us consider the onset of nematic order within the superconductor characterized by a gap 
amplitude A, and denote the critical value of s by s c (|A|). The essence of the picture of 
competing orders is that the onset of superconductivity should shrink the region of SDW 
order, and hence s° > s c (|A|). 

Then, we evaluate the critical point shift (s° — s c (|A|)) from coupling to the fermionic 
degrees of freedoms. It is determined by fermion susceptibility of the nematic order with 
and without SC. It is easy to obtain the susceptibility with SC, 



Xnem = / l(kf i A 10^/0 • ( 16 ) 
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(e(ky + \ 


A| 2 ) 3 / 2 



This integration looks proportional to magnitude of the SC order parameter, |A|. However 
the integral is non-zero as A — > 0, because of the finite density of states at the Fermi level. 

Here we expand the anisotropy term around the Fermi surface, and take the lowest term, 
('j(k) ~ |cos(2#fc)) for simplicity. Higher term corrections are subdominant, once we take 
the BCS limit(A/ J E F < 1). Then we obtain 

Xnem = ^^(l-a 2 |A| 2 ) , s c - s c (|A|) = — ^-^-a 2 |A| (17) 

where A/"(0) is a density of states at the Fermi energy, and a 2 is a cutoff-dependent positive 
constant. Therefore, the presence of SC suppresses the nematic order susceptibility, which 
indicates competition between the two order parameters. The a 2 term corresponds to the 
A term in the GL theory in Eqn. [TJ This calculation shows that a non-analytic term 
does not appear in the nematic order QCP shift , and the competition effect is relatively 
weak compared to the linear QCP shift as in spin density wave(SDW). 15 In other words, 




FIG. 1. Feynman diagram for the anisotropic coherence length vertex. The wavy line is for the 
nematic order parameter and the dashed lines are for the SC. The plain lines of the loop is for 
electrons's propagators with pairing. Note that directions of the fermion propagator is determined 
by the normal component of propagator. 

the competition basically follows the GL theory, so it is not necessary to keep fermions to 
understand the QCP shift of the nematic order parameter under SC. 

In the next section, we show that a non-trivial coupling appears in a finite momentum 
coupling corresponding to anisotropic SC coherence lengths by Fermi surface contributions, 
and the coupling is much stronger than usual GL couplings in the low momentum limit. 

IV. ENHANCED ANISOTROPY OF SC COHERENCE LENGTHS 

To incorporate the fermion contribution to anisotropic SC coherence lengths, we need 
to evaluate one Feynman diagram for the vertex term(r^*^) as in Fig. [TJ The electron 
propagator contains anomalous terms from SC pairing. The uniform SC order parameter 
is represented by |A|, and we consider the non-zero momentum component of SC modes. 
Also, in this paper, we only consider uniform nematic order parameter. In this case, it is 
easy to show that the vertex evaluation is equal to derivatives of the pairing susceptibility. 

Without loss of generality, we can assume that SC order parameter can be adjusted as a 
real number by choosing a proper gauge. In the SC phase, the pairing susceptibility depends 
on its direction. 




(18) 
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We define relations \xx = Xo — &X an d Xyy — Xo + $X so that 

pairing = J ^2 + + ^ + . . . (lg) 

If we assume that the band structure e(k) contains the nematic order parameter contri- 
bution, then it is easy to obtain the SC coherence anisotropy vertex by differentiating the 
pairing susceptibilies with nematic order parameter. Assuming e(k) is isotropic, the nematic 
order parameter just renormalizes a dispersion relation as follows. 

e(k) ->• e{k) + (p-?(k 2 ) cos(29 k ), (20) 

where ^{k 2 ) is for an isotropic magnitude dependent part. Then, 

pairing = r ^^|2 + f^(^ + ^*) + . . . ( 2 1) 

Full expressions for the vertex functions are 

r - dx °(rM -of ~ l n -i £k +p £k \( £k dsk \ i A o ( dek €k +p \ 

d<P [m=0 Z J k (E k+P + E k y [i + E k+p E k )[ E k d<P j + E k+P + E k [ d<P E^J 

£ = ^ X (p)|^ u = 2 / - 1 _ A o (^ggfc) + - A o ( ^ ek y (22) 
•90 A (-Efc+p + -Efc) 2 E k+p E k E k d(f) E k+P + E fc 90 E k+p El 

Here, we only focus on the BCS limit (Aq/ Ep — > 0). 

So far, our setup is very general. But to extract more information, let us introduce three 
parameters (Ci,C 2 , C 3 ) for the electronic density of states, effective mass, and the nematic 
coupling. With a general band structure, we can assume that momentum shifted dispersion 
relation is 

p2 

e k+p = e k + v F p cos(9 k - 6 P ) + C 2 — , (23) 

where C 2 is a non-universal number determining the effective mass. In general, the Fermi 
velocity has angular dependence with the four-fold rotational symmetry, and nematic order 
parameter breaks the symmetry down to the two-fold one. In this paper, we assume the 
Fermi velocity is Vp is isotropic and set to unity. Its angular dependence only changes 



numerical factors of the final results, specifically the constants (Ci, C 2 , C 3 ) of Eqn. 29 
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Also, we assume that the density of states at the Fermi surface is constant 



N{E F ) = dE F , (24) 

where C\ is a band structure related constant. 
Then, the above vertex functions are 

r L\_9rp f— f ^ ~ l n i £k +p £k \r £k d£k \ , A o i dEk €k +p \ 
WW ~ J 2n J de {Ek+p + Ek)2 U + Ek+pEk ^ Ek d(f) ) + Ek+p + E] }d<j> E k+p El> 

r ( ' dd [ ' A r ~ l A o ( £ kd£ k -Ag ggfc e fc 

1 F 7 2^7 d£ (£ fc+P + £,o 2 iW^ l £^^ ^ 

The energy integration requires a cut-off, A, and we assume it is much bigger than SC 
pairing. For the nematic coupling, we write 

-J^ = gi (k 2 )cos(26 k ), (26) 

where our third parameter appears in 7(/c 2 ) 

l^) = l + C^. (27) 
hip 

For fermions with a circular Fermi surface (e(k) = k 2 /(2m) — Ep), the introduced constants 
are 

C x =- C 2 = \ C 3 = l. (28) 

7T 4 

At the BCS limit, the gap function becomes negligible, but it is dangerous to set it to be zero 
in the integration because quasiparticles' pole structures are affected by the gap magnitude. 

Note that we do not place any constraints on p/A. By evaluating the vertex functions 
at the first oder in A/Ep, we find that 

r^(p) = 2C l p ^^{c 2 g 1 {p/A) + c 3 g 2 ( P /A)) 

2 2 

fV^(p) = 2cfi^^(C 2 g z {p/A) + C 3 S 4 (p/A)). (29) 

It is not easy to get analytical forms for the vertex functions, so we illustrate the Qi functions 
in Fig. [2j By performing asymptotic expansions, we obtain two limiting cases. The first 
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Si(x) g 2 (x) 




FIG. 2. Dimensionless functions of the vertex functions. See Eqns. 29 
limit is a metallic limit, p/A — > oo. 

hrM -> (3°) 

The other limit is a deep SC limit, p/A — )■ 0. 

^H>*M ~P 2 /A 2 

t^(p)~p 2 /A 2 . (31) 

Several remarks are in order. First of all, the Fermi energy scale from density of states in 
the vertex functions dropped out as expected. This is because we are considering a non-zero 
angular momentum channel in the vertex function, so the zero-momentum channel cancelled 
out. Next, the introduced three constants depend on microscopic details significantly. For 
example, different band structures from different lattices(or hopping parameters) change 
Fermi velocity and density of states, but one can absorb such differences into redefinition 
of the constants. Moreover, different order parameter structures such as different SC pair- 
ings(or different origin of nematicity) will affect the T^*^ function evaluation. However, 
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it is expected that final results will have the same functional behaviors with the adjusted 
constants. In other words, the functional behavior does not depend on details of mecha- 
nisms; different pairings like d-wave or s-wave and origin of nematicity only changes the 
three constants. 

In Fig. |3j we illustrate the vertex functions for the spherical Fermi surface case. 

p 2 — p 2 

r^*v(p) = CtFoip) — - 

f ^(p) = dJi(p)^^ (32) 

As we can see, the functional behavior of the anisotropy vertex is almost monotonic. The 
non-monotomic contribution comes from function which is much smaller than g\. From 
that, we understand that decreasing SC indicates increasing anisotropy effect. It certainly 
shows competition between the nematic order and SC. 

Finally, let us estimate how much the non-analytic term makes the anisotropy enhanced 
inside the vortex core compared to the usual GL term contribution. For the sake of simplicity, 
let us assume SC is completely suppressed inside the vortex core, and has full gap size outside 
the vortex core. Of course, its size is determined by the coherence length (gap size). Inside 
of the core, the functional behavior is approximated as 

Finside ~ \(pl -P 2 y ) (33) 
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Then, the enhanced anisotropy inside the vortex core, rj, is basically 

V = f d 2 p± ~ log(|) (34) 

up to some constant factors. We use a lattice length scale for comparison. This simple 
analysis shows that interplay between the Fermi surface and nematic order enhances the 
anisotropy of superconducting coherence lengths. 

V. DISCUSSION AND CONCLUSION 

Before closing the paper, let us take one step back and think about the meaning of our 
calculation. In theoretical physics, one of the most powerful concept is a low energy effective 
theory. It is very powerful because higher energy degrees of freedom are ignored and the 
theory become significantly simplified. Following the spirit, in a s wave SC, it is common to 
forget about fermion excitations in low energy limits because of the pairing gap in fermion 
spectrum. So, in our system, it seems reasonable to use a GL theory to describe systems 
because fermions are gapped. 

However, our calculation explicitly shows that such consideration is not enough and it is 
dangerous to use the naive GL theory. How can we understand such inconsistency? The 
answer is hidden in Fermi surface physics especially when multiple order parameters are 
considered. Even though the pairing makes fermion spectrum gapped, the SC order and the 
nematic order are coupled to each other with a finite momentum transfer. Once the finite 
momentum is bigger than the gap magnitude, fermions start to see Fermi surface physics 
and the energy gap becomes irrelevant. In terms of our calculation, it exactly corresponds 
to the metallic case, p£ 3> 1. 

To see importance of Fermi surface physics further, let us consider SC pairing suscepti- 
bility with nematic order parameter in the metallic ground state. 

X*(2p) = I , 1 , (9(-et +p ) - 9(et_ p )l (35) 

where the superscript (0) indicates presence of the nematic order parameter in band struc- 
tures. The denominator basically describes particle-particle channel energy difference and 
the step functions are for allowed phase space contribution. If we take p — > limit, then 
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the usual BCS logarithmic instability shows up, but the presence of the finite momen- 
tum regularizes the susceptibility. The anisotropy SC coherence length vertex is basically 
r</>V</<* = ^X^ / 9(p\(/,=o, and it is easy to check it indeed corresponds to the form of the metallic 
limit. 

Based on these considerations, the logic that gapped fermion theory only gives a GL 
theory is not correct. One should keep in mind that possibility of Fermi surface effects, 
especially in the case where a momentum (or frequency) dependent function is present. Such 
non GL type term is not universal, and depends on order parameter properties significantly. 

It is interesting to compare our results with recent numerical works on the Hubbard 
mo del was reported that even small anisotropic hopping term with strong interaction 

induces strong enhancement of anisotropy in some physical quantities such as conductance. 
The enhancement basically comes from interplay physics between SC, anisotropic hopping, 
and fermion excitations, and it seems the numerical calculations are consistent with our 
results. It would be very interesting if one can find a mapping between the current work 
and reported numerical results. 

So far, we mainly focus on a low energy theory to understand general competition effects 
between the two order parameters. Let us briefly comment on microscopic mechanisms on 
origin of nematicity in FeSe. First, it is unlikely for structural orthorhombicity alone to be 
an origin of strong vortex shape nematicity because it is too small even with non-analytic 
contribution from Fermi surfaces. Thus, electronic nematicity is likely to happen. Possible 
origins for the electronic nematicicty were suggested in literatures] 1 * 17 * 22 ! Magnetic ordering 
could be a candidate, but it was reported that there is no hole pocket and the nesting wave 
vectors ,(vr,0) or (0,7r), are supressed in FeSe^^ Alternatively, orbital ordering does no 
require any nesting, so it was proposed as another candidate for the nematicity. However, 
it is still an open question and remains to be understood in future works. No matter what 
is the microscopic mechanism of the nematicity, our main point is that magnitude of the 
nematic order is not required to be very big to explain the strongly distorted vortex shape. 

In this paper, we study interplay physics between nematic order and superconductivity. 
We show that competition between the order parameters is well captured by incorporating 
Fermi surface excitations coupling to the order parameters. It is shown that anisotropy of 
superconducting coherence lengths is descriebd by the non-local term induced by fermion 
excitations. We show that the non-local term becomes more important in metallic limit, 
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which shows (p^—Py) /p 2 behavior in the SC coherence anisotropy vertex function. Therefore, 
we argue the anisotropy channel becomes more enhanced near the vortex core, where SC is 
significantly suppressed. 
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